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ABSTRACT 

Dynamical friction arises from the interaction of a perturbcr and the gravitational wake it 
excites in the ambient medium. We study the effects of the presence of a boundary on dynamical 
friction by studying analytically the interaction of perturber with uniform rectilinear motion in a 
uniform homogeneous medium with a reflecting planar boundary. Wake reflection at a medium's 
boundary may occur at the edges of truncated disks perturbed by planetary or stellar companions 
as well as in numerical simulations of planet-disk interaction with no-outflow boundary conditions. 
In this paper, we show that the presence of the boundary modifies the behaviour of dynamical 
friction significantly. We find that perturbers are invariably pushed away from the boundary 
and reach a terminal subsonic velocity near Mach 0.37 regardless of initial velocity. Dynamical 
friction may even be reversed for Mach numbers less than 0.37 thereby accelerating instead of 
decelerating the perturber. Perturbers moving parallel to the boundary feel additional friction 
orthogonal to the direction of motion that is much stronger than the standard friction along the 
direction of motion. These results indicate that the common use of the standard Chandrasckhar 
formula as a short hand estimate of dynamical friction may be inadequate as observed in various 
numerical simulations. 



Introduction 



Dynamical friction, the reaction force a perturber feels from the gravitational wake it excites in the 
ambient medium, is one of the fundamental processes in the formation of astronomical struc tures. Applica- 
tions of dynamical friction to coll isionless systems include mass segregation in star clust ers (Kim fc Morris 



2003 



Portegies Zwart et al. 2004 ). galaxy distribution in hierarchical cluster formation ( Frenk et al. 19961: 



Bullock fc Johnston 2005 ) , the decay of satellites in galactic halos and galaxy mergers ( Binnev fc Tremainel 
2008h (and references therein). Dynamical f riction also occurs in collisional particle systems and was ap plied 



1985 



Stewart fc Wetherilllll988h. the 



to the study of the growth of planetesimals (iHornung. Pellat. fc BargdL 

eccen tricity excitation of planetary embryo orbits ( Ida 1990l:lKokubo fc Idalll99(f Namouni. Luciani. fc Pellat 
1996 ). and the confinement of planetary rings ( Goldreich fc Tremainel 19821) . For gaseo us systems, a 



Faltenbacher et al 



2005 



plications include the orbital decay of compact stars orbiting supe rmassive black holes (INaravanl 200^, 
the heating of intercluster gas by decelerating supersonic gala xies ( El-Zant. Kim, fc Kamionkowskil 2004 : 



Kim. El-Zant. fc Kamionkow ski 2005 ). the characterization of the differences in 



X-ray emi ssions in galaxy cluste rs between the standard cold dark matter model and modified newtonian 
dynamics (j Sanchez- Salcedoll2009h . 



Boundaries appear in dynamical friction because both the perturber and the ambient medium that 
surrounds it have finite sizes. In the common analyses of dynamical friction, these two boundaries enter the 
force through the Coulomb logarithm that requires cutoffs at small and large distances from the perturber 
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to avoid divergence. In the context of collisionl ess systems, the class ical derivation of dynamical friction 
using a pcrturber in uniform rectilinear motion (|Chandrasekharlll943f ) leaves the choice of the two cutoffs 
somewhat indeterminate and gives the force as: 



-167T 2 (GM) 2 TOl0g 



du u 2 f(u) 



V_ 

W 



(1) 



where G is the gravitational constant, M is the perturber's mass, m is the individual mass of a background 



medium star, V is the perturber's velocity, / is the medium's velocity distribution function, 



and 



are 



the small and large distance cutoffs respectively. The small distance cutoff, r m j n , is believed to be the larger 
of the perturber's radius and the distance related to the typical velocity in th e ambient medium wherea s the 
larger cutoff radius, r max , is usually taken as the size of the system (e.g. see (jBinnev fe Tremamdl2008l )). A 



more ac curate analy sis of the gravitational scattering of a star in rectilinear motion in a collisionless infinite 
system ( Henon|[l958h . shows that there is a natural large distance cutoff that equals the distance travelled by 
the perturber, Vt, where t denotes the time the perturber takes crossing the medium. This cutoff arises for 
durations smaller than the crossing time because the far stars do not have enough time to completely deflect 
the perturber by the standard (asymptotic) angle that depends only on the distance to the perturber and 
its velocity but not on time. The farther the star, the more incomplete the deflection, leading to a scattering 
angle that depends on time. Equivalently, it can be shown that only stars with impact parameters smaller 
than Vt substantially affect the perturber's ve locity. This analysis was a pplied to accurately estimate the 
relaxation time in an infinite dilute star system ( Ostriker fc Davidsenll 19681 ). When steady state is considered, 
the distance travelled by the perturber becomes comparable to the system's size and the upper cutoff, r ma x, 
is taken as the system's size. Similar prescriptions are applied to circular motion in spherically symmetric 
systems where it is shown that only the vicinity of resonances within the system contribute to dynamical 
friction (jTremaine fc Weinberg||l984l ) including the standard long and short distance cutoffs. 



For a perturber moving in a gaseous medium, the small distance cutoff is taken as the perturber's 
accretion radius or physical size. The l arge distance cutoff depends on the perturbat i on regime. In steady 
state , it is taken as the system's size (|Dokuchaevl 11964 iRuderman fc Spiegell Il971t iRephaeli fc Salpeter 
19801 ) whereas in the time-dependent regime, in which the perturbation's duration is much smaller than the 
crossing time, it is given by the distance travelled by the sonic Shockwave that the perturber excites i n the 
medium JOstrikerlll999l ISanchez-Salcedo fc Brandenburg|[l999[ l200ll : IrCim fc Kimll2007t IWolfsonl[l977t ) and 
( Namouni 2010h hereafter Paper I. Although the nature of the large distant cutoffs used to truncate the 
density perturbation is of different physical origin for collisionles s and collisional media, the corresponding 
friction forces are similar (ISanchez-Salcedo fc BrandenbursJl200lh . 



In both collisionless and collisional media, dynamical friction was determined in the two limiting cases 
of steady state and the time-dependent regime. For the latter, the independence of dynamical friction of the 
medium's size is straightforward as the perturbed region that acts back on the perturber is much smaller than 
the ambient medium. The situation is much less clear for steady state in which the density enhancement 
raised by the perturber has reached the medium's boundary. Suppose for instance that a perturber is 
set in motion near the boundary of a gaseous medium such that the distance the sonic Shockwave has to 
travel to reach the boundary is much smaller than the medium's size. Before the Shockwave reaches the 
boundary, dynamical friction is in the time-dependent regime. However as the boundary is reached, the 
medium's full size cannot serve as a large distance cutoff. Instead, the proximity of the perturbed region to 
the medium's boundary will affect the way the density enhancement acts back on the perturber depending 
on how the boundary responds to an incoming Shockwave. It is the object of this paper to find out how 
dynamical friction may be changed as the perturber's Shockwave reaches the medium's boundary. As such 
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a task is quite difficult, wc will make a number of simplifying assumptions aimed at illustrating that the 
classical behavior of dynamical friction may be changed drastically because of the presence of a boundary. 
We assume that the perturber moves in a homogeneous gaseous medium along a rectilinear trajectory with 
uniform velocity near the medium's boundary that is taken to be planar and reflecting. The boundary's 
curvature is neglected as the medium's size is supposed to be large compared to the perturber's distance to 
the boundary. Wake reflection at a m edium's boundary may occur at the edges o f truncated d i sks perturbed 



by p lanetary or stellar companions |Hahnl |2003[ Larwood fc Papaloizoul 



1997 



War. 



d 



2003; 



Tanaka et al 



2002 ) as well as in numerical simulations of planet-disk interaction with no -outflow boundary conditions 
(jCresswell et al.ll2007t [Nelson et al.lboool: IPapalo izou. Nelson, fc Massetll200ll) . Interestingly, we shall show 
that most of the significant changes on dynamical friction come not from the reflected component but from 
the truncation of the Shockwave incident on the boundary. To derive the friction force for this set-up, we 
shall use the general expression of dynamical friction derived in Paper I that allowed us to study the influence 
of a perturber's acceleration on dynamical friction. 

The paper is organized as follows: in section 2, we recall the general expression of dynamical friction 
derived in Paper I and characterize how the medium's size affects dynamical friction in general. In section 3, 
we derive the dynamical friction force for a perturber moving normal to the medium's boundary and show 
that perturbers are pushed away from the boundary and reach a finite terminal velocity. In section 4, we 
examine dynamical friction for motion parallel to the boundary and show the appearance of an additional 
force component normal the trajectory that pushes the perturber away from the boundary. Section 5 includes 
a summary of main results and possible further refinements. In Appendix A, wc characterize the gas flow 
forced by the perturbers in the configurations studied in sections 3 and 4 by deriving the corresponding 
velocity field and density perturbation and linking them to the known results of the Bondi-Hoyle-Lyttleton 
accretion. 



2. General dependence of dynamical friction on the medium's size 

The reaction force exerted on a perturber, with a general potential (f> p and a general trajectory £(t), by 
its gravitational wake in a homogeneous gaseous medium, of unperturbed density po and sound speed c, was 
derived in Paper I and reads: 

F = j<mGpof d3x g_ Pp{x _ Y+A)VMn (2) 

c J dV[X+$(t)],\Y\<ct l* I 

where p v is the perturber's density distribution, p p = V 2 p /47rG, A = £(t) — £(t — \Y\/c) and J$? (t) is the 
Heaviside function. The spacial integration domain dV(x) is defined by the medium's boundary as well as 
the lower distance cutoff in the perturber's vicinity such as its size or its corresponding accretion radius. The 
boundaries are used to truncate the perturber's potential at small and large distances to avoid the force's 
divergence. For a point-like perturber of mass M, the expression is simplified to: 

F= ^ / rsmUdrdvdip -—r^, (3) 

c J dV[y+£{t-r/c)\,r<ct \V ~ ^\ 

where r, 9 and tp are the spherical coordinates of the vector y defined as the relative position within the 
medium with respect to the perturber's position at the retarded time t — r/c. It may be written as y = 
x — £ (t — r/c) where x is the position of fluid clement in the medium. This definition of y influences the force 
integral through the boundary conditions of the medium denoted by dV. The phase A = £(t) — £(t — r/c), 



-4- 



thc upper radial cutoff at the retarded sonic Shockwave r < ct, and the boundary conditions are all that is 
necessary to compute the friction force. In particular, as the force expression is independent of the density 
perturbation, no explicit expression for the density enhancement is required for the force calculation (for 
more details see Paper I). The derivation of expression ([3]) assumes that the perturber is set in motion at 
t = a fact that introduces the size of the retarded sonic Shockwave as an upper cutoff. This assumption is 
not restrictive as it allows us to study the time-dependent evolution of dynamical friction as well approximate 
steady state behaviour. It models a perturber that just entered the ambient medium such as a star cluster or 
a dwarf galaxy that fall into a larger galaxy. It also de scribes turning on a numerical simulation of dynamical 
friction ( Sanchez-Salcedo fc Brandenburg 19991 2001 ). 



The formula ([3]) allows us to see more clearly how the classical formula of steady state dynamical 
friction arises and how its dependence on the medium's size is erroneous. Consider a perturber with uniform 
rectilinear motion that travels in a medium of a typical size r max . Writing the trajectory as £(t) = Vte z 
where e z is the unit vector along z, the phase A = ^#re 2 where jtft is the Mach number. The friction force 
is equally along e z and its expression is simplified as follows: 



F 
F 



Jf(t)(GM) 2 Po 



c JdV[y+£(t-r/c)],r<ct 

2irJ?(t)(GM) 2 p r dr 

V 2 JdV[y+£{t-r/c)],r<ct r 



r sin 9dr d9 dip 



• cos 9 — A 



r 2 + A 2 -2rAcos6>| 3 / 2; 



(1 + 



-y/2 



(4) 
(5) 



where in the last expression the integrals over p> and r 
and T max are determined from the boundary conditions. 



cos 9 were performed. The extremal values r„ 



The condition on the medium's extent may be written as \x — < r max defining a sphere that comoves 
with the perturber. Lastly, we supplement the integration domain with a lower cutoff ro, usually taken to 
be the perturber's size or accretion radius, 2GM/V 2 , and implemented as \x — £(t)\ > r . Expressing these 
conditions in terms of the variables r and r that appear in equation ([5]) shows that we need to find the 
domains that satisfy: (i) -1 < r < 1, (ii) r < r(l + Jt 2 - 2.^r) 1 / 2 , (iii) r(l + ./# 2 - 2./#r) 1 / 2 < r max , and 
(iv) r < ct. The first three conditions yield the following integration domains: 



ro 



+ 



< r < 



< r < 



< r < 



ro 



with r min = -1, T max = 



1 + Ji 2 - (r /r) 5 



2Jt 



With T min = -1, r max = 1, 



1 



1 + ^ 2 - (r raax /r) 2 

With T mirl = ——, , T„ 



= 1. 



(6) 
(7) 

(8) 



The fourth condition depends on the distance ct travelled by the retarded wavefront; it sets the dynamical 
friction regime: time-dependent or steady state. In the time dependent- regime, (ct <C r max /|l + ^\), 
the sonic Shockwave has not reached the medium's boundary whereas in the steady state regime, (ct 3> 
?"max/| 1 — •^\)^ the peturbation has travelled the size of the medium. Calling / the integral that appears in 
([5]), the integration for the time-dependent regime extends over the two domains (|6|7j) with ct substituted 



for 



II 



and I is given as: 



log 



1 + . 



- [1 - sign(l - Jt)\ log 



f- |l-^T|][l + sign(l- 
(1 - M)ct 



'0 



(9) 
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(see Paper I fo r a detailed calculation) . This is the classical friction expression in the time-dependent regime 



(Ostrikerll999tlSanchez-Salcedo fc Brandenbur jl999ll200ll;lKim fe Kimll2007tlKim. Kim, fe Sanchez-Salcedol 



20081 ). For the steady state regime, |c — V\t is much larger than the size of the system. This condition ignores 
how the medium's boundary reacts to the incident Shockwave and how the latter is truncated, reflected or 
transmitted by the boundary. As ct > r max /|l — the condition from the propagation of the retarded 
wavefront (iv) becomes irrelevant and the three domains (|6I7I8|) contribute finite terms to the force. The 
inner and middle domains yield the same result as in (j9j) except that ct is replaced with r max /|l + The 
outer integration domain ([8]) yields the term: 



1 + M 



1 



- i [1 + Ji - |1 - Jf\] [1 + sgn(l - J?)] 



which, when summed up with the terms of the first two domains, yields exactly: 



sgn(l 



log 



(10) 



(11) 



By compensating for the contribution of the inner domain ([6]) , the force term Tjgj cancels dynamical friction 
for subsonic perturbers and yields a force on supersonic perturbers similar to the Chan drasekhar formula 
for collisionless systems. This is the classical result of steady state dynamical friction (jDokuchaevi 11964 ; 
Ruderman fe Spiegel 1971 ; Rephaeli fe Salpeter 198C)I) . The cancellation for subsonic perturbers happens 
only because of the choice of an outer boundary that is symmetric with respect to the perturber and 
comoves with it, i.e. \x — £(t)\ < r max ; the two conditions are necessary. The change of behaviour from the 
time-de pendent regime to steady state may be understood through the explicit expression of the perturbed 
density (jQstrikerlll999l : I Just fc Kegellll990h : 



g(x,t) 



GMk/c 2 



[(Z-Vt) 2 + R 2 (l-^ 2 )] 



1/2 



where k 



1 if R 2 + z 1 < (ct) 2 , 

2 if Jt> 1, R 2 + z 2 > (ct) 2 , 

Y2 1\V2 







z + R[M A - 1) ' < Vt, 
and z^i < ct, 
otherwise. 



(12) 



where R = rsin# is the cylindrical radius. Whereas the density perturbation is symmetric with respect to 
the perturber, the sonic wavefront of radius ct that was launched at t = is not symmetric with respect 
to the perturber as it does not comove with it (Figure 1, top left panel). In the intermediate regime, when 
ct <C r max , the asymmetric part of the density perturbation inside the sonic wavefront results in dynamical 
friction. As the wavefront reaches the comoving symmetric boundary, friction decreases and tends to zero 



when ct 3> 



as the density perturbation artificially acquires front-back symmetry with respect to the 



perturber (Figure 1, bottom left panel). However, this explanation that is based on the specifics of the 
density perturbation works only when uniform motion is considered as in general, motion is accelerated (or 
decelerated), as by dynamical friction itself, and the density enhancement generated by accelerated motion 
nowhere has front-back symmetry (Paper I). Discarding the unphysical assumption of truncating the force 
at large distances using a comoving symmetric boundary (which is equivalent to implicitly assuming that 
perturber is at the center of a symmetric medium) will lead to significant changes in dynamical friction as 
suggested by Figure 1 (right column) where the surfaces of equal density are shown for a subsonic perturber 
in a medium with an arbitrary shaped boundary that does not comove with the perturber. As the sonic 
wavefront reaches the outer boundary, the density enhancement becomes asymmetric with respect to the 
perturber and friction is not only finite but will have a component perpendicular to the perturber's trajectory. 
In the following two sections, we examine a simpler boundary and show quantitatively that the properties of 
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dynamical friction change significantly when the assumption of a comoving symmetric boundary is discarded. 
We also remark that most of our conclusions regarding the use of a comoving symmetric boundary apply to 
collisionless systems. In effect, for rectilinear motion, the derivation of dynamical friction from the velocity 
impulse caused by a test star on the perturber is invariably centered around the latter. When steady state 
is considered, the cutoff r max is used in the perturber-centered frame which is equivalent to assuming a 
comoving symmetric boundary. 



3. Force for motion normal to the boundary 

As discussed in the previous section and shown in Figure 1, dynamical friction in the presence of a 
boundary unrelated to the perturber has parallel and normal components to the perturber's motion. To 
study each component separately, we choose perturbers that move normal to the boundary (this section) 
and parallel to the boundary (section 4). We therefore consider an infinite medium with a reflecting planar 
boundary that occupies the space z > and a perturber of trajectory £(t) = (Vt + d)e z where e z is the unit 
vector normal to the boundary. This setup models dynamical friction on a perturber whose outgoing wake 
has reached the nearest part of the medium's boundary or one that has just entered a gaseous medium. The 
assumption of a planar boundary is valid if the perturber is not near the physical center of the medium. No 
other forces apply to the perturber in its motion in the medium. As the boundary is reflecting, the perturber's 
density enhancement will have two components: the outgoing wake incident on the boundary and truncated 
by it, and a reflected wake. The reflection assumption serves two purposes, first as the large distance cutoff 
(akin to the small distance cutoff) is commonly used to truncate the density enhancement and estimate 
dynamical friction, the outgoing wake truncated at the boundary will not be affected by reflection; the 
corresponding force would be the same for a non-reflecting boundary. Second, at little extra computational 
cost, we obtain the force caused by the reflected wake. Wake reflection at a mediu m's boundary occurs for 



1997 


Wardll2003: 


Cresswell et al.ll2007) as well as in numerical simulations of planet-disk interaction with no 


outflow boundary conditions ( Cresswell et al. 


2007; 


Nelson et al. 


2000[ 


PaDaloizou. Nelson. & Masset 


2001) 



The sonic Shockwave excited by the perturber reaches the boundary at t = d/c. Prior to this time 
(0 < t < d/c), the reaction force from the wake on the perturber is that of the time-dependent regime (0). For 
later times, dynamical friction remains along the direction of the perturber's motion because the medium's 
boundary has a rotational symmetry with respect to the perturber's trajectory. We consider separately the 
forces from the outgoing and reflected wakes and assume that the perturber is not close to the boundary, i.e. 
d (1 + ^#)ro/|l — ^#|. Analytic expressions as well as portraits of the density perturbation and velocity 
fields are given in Appendix A. For the outgoing wake of a perturber moving away from the boundary, the 
geometrical truncation of the density perturbation will reduce dynamical friction. As the surfaces of equal 
density arc symmetric with respect to the perturber, truncating part of the wake on the side opposite the 
perturber's motion (Appendix A, Figure 9) will effectively reduce the drag behind the perturber. In fact, 
we show below that the situation is more dramatic for slow subsonic perturbers as truncation reverses the 
direction of the force. The reflected wake remains behind the perturber for supersonic motion and overtakes 
the perturber at t = 2d/(l — ^) for subsonic perturbers. Before this time and at all times for supersonic 
motion, the reflected wake is simply an additional density enhancement behind the perturber and increases 
the drag above that felt in an infinite medium because the same part of the wake is now closer to the 
perturber. After the reflected wake has overtaken a subsonic perturber, the amplitude of its drag decreases. 
We show below that its friction force may be reversed for slow subsonic motion in manner similar to that 
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of the outgoing wake's force. Dynamical friction is derived for times d/c <C t <C r max /c. The first condition 
defines a local steady state with respect to the boundary nearest to the perturber and the second condition 
ensures that we do not follow the perturber's wakes (outgoing and reflected) as they reach and bounce off the 
whole (closed) physical boundary of the medium. For a perturber moving towards the boundary, dynamical 
friction is derived for times t < d/V prior to the perturber's exit from the medium. Only subsonic motion 
is considered as the boundary ahead of a supersonic perturber does not affect dynamical friction. 



3.1. Force from outgoing wake 



Wc first assume that the perturber moves away from the boundary (i.e. > 0). The boundary 
condition z > written in terms of the variables of expression ([5]| becomes r > j$ — £(i)/r with — 1 < r < 1. 
The condition that the volume that contributes to the force be that inside the retarded sonic wavefront, i.e. 
r < ct, then implies that T min — £(i)/(l +^#) and r max = 1 only for t > d/c. The integration domain related 
to the lower cutoff radius, r , ([5]) remains unchanged. For larger radii, we have: 



ro 



at) 



< 



< 



< r < ct with 



with T min = -f, 

_ _ , m 



= i, 



t. 



(13) 
(14) 



The domains (|6]) and (fT3|) contribute force terms identical to those of equation (|9|) except that ct is replaced 
with £(i)/(l + ./#). The third domain (fT4|) related to the truncation of the sonic wavefront contributes the 
new term: 

^sgn(l - JK) (l-.J? 2 )r + 



7 JH 



i ■ // 



dr. 



(15) 



r r 3/2 [(j _ 

Inspection of the boundary conditions associated with a perturber that moves towards the boundary < 0) 
shows that the term I rj^ equally applies to — 1 < < provided that the perturber does not leave the 
medium (t < (d — r )/|F|). Summing up the contributions of the three domains (|6I13I14[) gives: 



I 



(1 + + IJtdjct) 2 + Jt - f + 



lOB 



^f 2 ) 2 log 
d/ct + J( 



(d/ct + J() 



i 



\ + 

for 



ct + [(I + J/ 2 + 2JZd/cb){l - Jt 2 )\- 
\Jt\ < 1, 



(16) 



I = -3 (I + .£ 2 + IJCdjct) 2 + ZJi + 3 + 



sm 



2Jt 



— sin 



- I 



2Jt[J{ + d/ct) 



los 



(Vt + d)(^Z - l)ct 



for Jt >\. 



(17) 



For the local steady state with respect to the boundary, t ^> d/c, the terms d/ct disappear yielding dynamical 
friction as a function of with the exception of the Coulomb logarithm for supersonic motion. The force 
described by expressions (|I6II7p is shown in Figure 2. As expected the truncation of the outgoing wake 
reduces the amplitude of the force for subsonic as well as supersonic perturbcrs that travel away from the 
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boundary. For subsonic perturbers moving toward the boundary, the truncation of the wake region ahead 
of the pcrturbcr increases the dynamical drag as more weight is given to the density enhancement behind 
the perturber. For —1 < ^# < 0.43, the truncation of the wake leads to the acceleration of the perturber 
away from the boundary. Motion at ,M = 0.43 is in stable equilibrium as for larger (smaller) velocities the 
force decelerates (accelerates) the perturber. Another interesting effect of the boundary's presence is the 
instability of perturbers at rest in the medium. While the integral / (jT6|) vanishes as tends to zero, the 
force which is proportional to I j M 2 does not. It is given as: 

fh » d/c, |ur| « i) = _^m^£^iMl. (18) 

Perturbers at rest are unstable to small perturbations (Figure 5, bottom left panel) and tend to be accelerated 
toward the terminal velocity V = 0.43c. This instability originates in the fast expansion of the truncated 
sonic wavefront with respect to the perturber's initial velocity (c/V 3> 1) giving most of the weight in the 
reaction force to the matter ahead of (behind) a perturber that moves away from (towards) the boundary. 



3.2. Force from reflected wake 

In Appendix A, it is shown that the density perturbation associated with the wake reflected by a 
planar boundary is generated by a fictitious pcrturbcr whose trajectory is the mirror image of the actual 
perturber's motion with respect to the boundary, x ■ n = 0. The fictitious perturber's trajectory is £ r (t) = 
£(t) — 2[£(t) • n}n, where n is the normal to the boundary. The force, however, is applied at the perturber's 
position, Its expression is therefore given as: 

,^(t)(GM) 2 p f 3 ^ 5[u-t+\x-C(u)\/c]jr(u) x-£{t) 



F c 2 J dv{x) dxdU \x-?(v)\ l*-£WI 3 ' 

F = f rainMrd^J^, (19) 

where y = x — £ r (t — r/c), r, 9, tp are its spherical coordinates and S = £(t) — £ r (t — r/c). For uniform 
rectilinear motion along n — e z , £ r (t) = —£(t) = —(Vt + d)e z , y = x + £(t), and X = [2£(t) — ^r]e z . The 
force is along e z and its amplitude is given as: 



F = 



2ir(GM) 2 J^(t)p Q f dr 



V 2 .lav r 



Jt 2 (\ - Sr/r) 



(E/r) 2 [l + (S/r) 2 -2Sr/r] 1 /2 



(20) 



The three boundary conditions, \x — £(t)\ < r , z > and r < ct show that the force is evaluated differently 
for subsonic and supersonic perturbers. There is only one integration domain as long as the wake does not 
overtake the perturber. This applies to supersonic perturbers at all times and only for t < (2d— ro)/(l — ^M) 
in the case of subsonic perturbers. This domain is given as: 

: < r <ct with r min = - r max = 1. (21) 

1 + M r 



1 lt may be questioned whether the terminology 'dynamical friction' applies to the force that accelerates perturbers instead 
of decelerating them. We choose to retain 'dynamical friction' as the accelerating force derives from the same phenomenon and 
equation as standard dynamical friction. It is only the boundary's presence that is responsible for the behaviour change. 
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Calling J the integral in equation ([20]) , these boundaries yield: 



J = 



dr 



i . * 



sgn(l-£/r) (1 - J( 2 )r 2 + Z£Jlr - 2£ 2 



r(E/r) 2 E 2 r 1 /2[(i_^2) r + 2 ^]i/2 



(22) 



J = - lo 



{l + 2d/Vt){JZ + \) 



(1 + d/Vt)(„¥? + 2) 
Jt(\ + .J( 2 + 2J?d/ct) 1 / 2 



2(1 -d/ct) J£(J( + 1) 

+ 2)(l + 2d/W) + J{ + 2 



los 



^ + 2d/cf 

.#[1 + (1 + + 2^/ci) 1 / 2 ] 



Jo 



Jo 



Jn = 2 



(^T + 2)[1 - (1 + ..# 2 + 2^d/ct) 1 / 2 ' 

1 + ^Trf/ci + [(1 - .^ 2 )(1 + J( 2 + 2^/ci)] 1 / 2 



(23) 



(l-^f 2 ) 1/2 lo 

_ 1)1/2 



+ d/ct) 1 + (1 



1/2 



, for J( < 1, 



sin 



- 1 



2J£{JC + d/ct] 



— sm 



~2~£ 



for > 1. 



In deriving these expressions we used the fact the first term in equation (|22p is negative as E > r for 
r < 2£(i) /(l + The wake completely overtakes the subsonic perturber after t = (2d + ro)/ (1 — ^#). In 
this case the boundary conditions yield three domains as follows: 



gg(t) - r 

\+..# 
2^(*) + T-o 



< 



< 



< 



2j(t) - r 

< — Wltn T„ 



< 



1 • .// 



with r n 
with r„ 



r 
r 

m 

r 



1, 



2rE 



= 1. 



(24) 
(25) 
(26) 



The middle domain (|25p evaluates the force near the perturber's location and it can be easily shown to 
contribute a force term that vanishes as ro tends to zero. Calling J\ the integral over the first and third 
domains (|24l26p defined by equation (f2T)|) . and noting that the sign of the first term in (12"2"|) is positive in the 
third domain, we find that the force integral on a subsonic perturber overtaken by the reflected wake is: 



Ji = J + 21og 



(J( +!)(! + 2d /Vt) 



2(1 + d/Vt) 



2Jt(\ - Jt - 2d/ct) 
J£ + 2d/cb ' 



(27) 



As for the outgoing wake, inspection of the boundary conditions for a perturber moving towards the boundary 
shows that the expressions (|23l24l27p are valid for < 1. The force from the reflected wake is shown 
in Figure 3. For supersonic and subsonic perturbers ahead of the wake, the reflected wake increases the 
dynamical drag. For larger Mach numbers, the contribution of the reflected wake is much smaller than 
that of the outgoing wake. The force amplitude is reduced after the reflected wake overtakes the subsonic 
perturber. Below ^# = 0.27, the reflected wake accelerates subsonic perturbers in the same manner as the 
outgoing wake. Perturbers at rest are equally unstable and the local steady state force from the reflected 
wake for slow subsonic perturbers is equal to that from the outgoing wake (|18D . 



The reaction force from both outgoing and reflected wakes is shown in Figure 4. Dynamical fiction is 
reversed at M = 0.37, a velocity that defines the effective stable rest state. The action of the reflected 
wake makes the dynamical drag beyond ^# = 0.57 stronger than that in an infinite medium peaking near 
M = 0.8 at 20% above that of equation ((9j> . Zero- velocity perturbers are unstable and are subjected to a 



- 10 - 



force twice that from the outgoing wake (|18[) . As a result of the boundary's presence, supersonic perturbers 
feel a smaller drag force that is about 10% below the value of the force in an infinite medium. The peak of 
dynamical drag remains near Mach 1. 



For a perturber moving in an infinite homogeneous medium occupying the space x > 0, parallel to its 
reflecting boundary, the gravitational wake will be truncated only on one side parallel to the perturber's 
trajectory £(t) = de x + Vte z (Appendix A, Figure 10). Motion parallel to the boundary is therefore subject 
to an additional component normal to the trajectory unlike motion normal to the boundary which results 
only in a force along the perturber's trajectory. As t becomes much larger than d/c, the forces from the 
outgoing and reflected wakes along the perturber's trajectory will tend to similar values close to half that of 
expression © as the wake is approximately cut by half (Figure 10). The total force along the trajectory will 
therefore tend to its value in an infinite medium. The force normal to the trajectory however will increase 
with the size of the perturbed region and force the perturber to move away from the boundary. 



The force from the outgoing wake is derived using the general expression ([3]) . With our choice of medium 
with x > and = de x + Vte z , A = ^re z . The main difference in the derivation of the force with 
respect to sections 3 and 5 is the integration domain of the variable ip that not longer runs from to 2ir as 
rotational symmetry is broken for r > d. Dividing the integration domain over r into Tq/\\-\-^\ <r<d and 
d < r < ct and writing F = F x e x + F z e z , the force related to the former domain F x (ro/\l + ^f\ < r < d) = 
whereas F z (r /\1 + < r < d) is given by equation §§§ where ct is replaced with d. In writing this we 
assumed that d 3> (1 + ^#)r /|l — so that the perturber is not too close to the boundary. For the part 
related to the domain d < r < ct, the force expressions may be written as: 



where the integration over ip was performed. In the variables of equation the boundary condition x > 
becomes rsin^cos^ > —d. Denoting, cosipd = — d/rsm9 and sinQd = d/r, the boundary condition yields 
the angles ip m in and ip max as functions of 9 and r as follows: 



4. Force for motion parallel to the boundary 



4.1. Force from outgoing wake 




(28) 



(29) 



^min 



0, <p : 



max — 



2tt for < 9 < 9 d 
ifd for 9 d <9 < 7r 



and 7r — 9d < 9 < ir, 



0, 



'd- 



Substituting these expressions into equations (|28I29|) yields: 




r2W2 



- 11 - 



F z {d <r<ct) 



Jf(t)(GM) 2 p 



2 + 2. // COS u j 
2 „ pet rn-6 



3/2 



2(cos0- Jt) sin0d0dr 



+ 



c " J r =dJe=e d r (1 + ^ 2 - 2^#cosi 

Jf(t)(GM) 2 p Q f ct f 0d 



dJ9=0 



,3/2 



COS? 



(1 + .-# 2 - 2^ cos ( 



^3/2 



COS 



(1 + . //-' + 2JV cos I 



^3/2 



(31) 



These expressions cannot be reduced analytically if local steady state is not assumed. To follow the force's 
time evolution as it reaches steady state, the expressions are integrated numerically, combined with the force 
term associated with the domain ro/|l + j4t\ < r < d and shown in Figure 5. As expected, it is apparent on 
the magnitude of F z for subsonic perturbers and on the slope of F z (t) for supersonic perturbers that the force 
along the trajectory F z tends to half the value given by the force on a perturber in a infinite medium (j9]). We 
note that while this result is straightforward when looking at Figure 10, it could not have been guessed from 
the formula used for an infinite medium (j9]) as only the constant unperturbed density represents the medium 
in that expression and it is not reduced by half if there is a boundary. The force normal to the direction 
of motion, F x , appears as soon as the boundary breaks the rotational symmetry of the wake incident on it 
after t = d/c. The force from the reflected wake always pushes the perturber away from the boundary and 
increases linearly with the size of the perturbed region which is proportional to log(ct). It is strongly peaked 
at ^ = 1 where it reaches its finite maximum. In the local steady state regime, ct 3> d, the normal force is 
given by the first term of equation (|30[) which may be simplified to: 



F,. 



2jr{t)(GM) 2 p \og(ct/d) 
V 2 (l + y£) 



(1 + Jt 2 ) K 



4^ 



1 



(1 +JZ) 2 E 



(1 + Jif 



(32) 



where K and E arc the complete elliptic integrals of the first and second kind respectively. This expression 
is not valid at = 1 as it gives an infinite force (the maximum amplitude of F x (|30j) at ./# = 1 grows 
faster than (log ct/d) 2 and its value becomes infinite in an expansion in terms of d/ct). For slow subsonic 
perturbers, ^# <C 1, the force is independent of velocity and may be written as: 



Fr 



TTjf{t){GM) 2 Po \og(ct/d) 



For fast supersonic perturbers, 



3> 1, the normal force is given as: 
nJ^(t)(GM) 2 cp \og(ct/d) 



F, 



(33) 



(34) 



This expression of F x bears some resemblance to that of F z © except that (i) the former's amplitude is 
smaller than the latter's by a factor (ii) it is the distance to the boundary that takes on the role of r Q in 
the Coulomb logarithm, (iii) the size of the perturbed region is ct instead of Vt for supersonic perturbers. 



4.2. Force from reflected wake 

The force from the reflected wake is derived using the general expression (flU)) with S = £(t)—£ r (t—r/c) = 
2de x + J%re. z . In the variables of equation (|19p , the condition x > becomes r sin cos Lp > d implying that 
d^TT0,i<6<TT — 9 c i and ipd — K^P^K — Pd- The boundary condition related to the lower cutoff distance 
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|y — S| > ?'o shows that there is a range in radius between r± ~ (2d sin 6 ± di)/(l + Jt 2 — 2Jt cos9) with 
d\ = rl{l + Jt 2 - 2 Jt cos 9) - Ad 2 (cos 9 - Jt) 2 where 

r 2 (l + Jt 2 - 2^cos6») + 4d 2 - rl , . 

cos <p < , . (35) 

4rd sin # 

This condition is the counterpart of the middle domain (|25j) associated with perturbers moving normal to 
the boundary that samples the force from the reflected wake in the perturber's vicinity. In the limit where 
ro tends to zero it can be shown easily that this condition does not affect the boundary values of the variable 
ip. The force from the reflected wake is therefore given as: 

J^(t)(GM) 2 p f ct r-K-Vy=*-va 2r(rsin0cos^-2d)sin0d0drd</> 
F x = ^ 



F- 



r =dJ8 d J v =o [r 2 (l + Jt 2 -2Jt cos9) + 4d 2 ~ 4rdcosipsm9} 3 / 2 

(36) 

J^(t)(GM) 2 p f ct r-x-Otr'p^-'Pd 2r 2 (cos 9 - Jt) sin 9d9drdip 

c 2 Jr=dJe d J<p=o [r 2 (l + Jt 2 - 2Jt cos 9) + 4d 2 - 4rd cos ip sin 9} 3 / 2 

(37) 

Like those of the previous section, these expressions cannot be reduced further analytically. They are 
integrated numerically to see how the force reaches local steady state and shown in Figure 6. The component 
along the trajectory, F z , tends to a value slightly smaller than half the value of ((9]) because of the position 
shift d with respect to the outgoing wake - for supersonic motion, the factor half is seen on the slope of 
the force. Local steady state is reached after a transitory period where the reflected wake is located in the 
vicinity of the boundary and does not cross the perturber's trajectory. This is the reason why F x initially 
pulls the perturber to the boundary before reaching its asymptotic positive value that is proportional to 
log ct. Like the force from the outgoing wake, F x is maximal at ^# = I. The analytical expressions of the 
normal force from the reflected wake in the local steady state regime, ct ^> d, are identical to those of the 
outgoing wake (|32I33I34[) . 

The reaction force from both outgoing and reflected wakes is shown in Figure 7. The force along the 
trajectory is comparable in amplitude to that in an infinite medium except in the neighborhood of Mach 1 
where its is larger by about 15%, a trend that is opposite to that of the friction force acting on a perturber 
moving normal to the boundary. The force normal to the trajectory increases linearly with \og(ct/d) (for 
Jt 7^ 1) and is maximal at Mach 1. For < 1, F x > \F Z \ whereas for Jt > 1, F x <C |-F*|. We note that 
like for motion normal to the boundary, the gravitational wakes tend to push the perturber away from the 
boundary. 



5. Conclusion 

We set out to examine the effect of the presence of a boundary on dynamical friction. We showed 
that the expressions of the force in infinite (j9]) and finite (fTTj) media may not be used in the presence of a 
boundary such as the medium's outer and inner limits or any other internal edge. In particular, dynamical 
friction along the direction of motion may be reversed and a new drag component perpendicular to the 
boundary is present and forces the perturber to move away from the medium's edge. These results are of 
particular relevance for the numerical simulation of dynamical friction once the perturber's wake reaches the 
boundaries of the physical system or the simulation box as for instance in problems related to disk-planet 
interactions and galaxy mer ging where timescales are ob served to differ significantly from the prediction of 



the Chandrasekhar formula (|Bovlan-Kolchin et al.ll2008l ). 
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We have considered separately motion normal and parallel to the boundary to isolate the effects of 
the instability of zero- velocity perturbers, the reversal of dynamical friction, and normal dynamical friction. 
Whereas qualitatively accurate for a perturber on a trajectory inclined by an arbitrary angle with respect 
to the boundary these effects may not be added quantitatively linearly because dynamical friction depends 
on the phase vectors A and X. This is why it is not straightforward to carry over the first two effects to 
rotational motion in bounded systems (such as disk-planet interactions) as they would require the pcrturber's 
orbit to have some eccentricity and hence a velocity that makes a time-dependent angle with the boundary. 
Normal dynamical friction however readily applies to circular motion in a spherical system. The dependence 
of the normal force on the distance to the boundary is translated for a perturber in circular motion around an 
attracting center as a radial force that depends on the orbital radius and the location of the medium's edge 
for subsonic perturbers, and additionally on the velocity for supersonic perturbers. This time-dependent 
force affects the pcrturber's energy but not its angular momentum and is most efficient where motion is 
resonant with the sound speed. 

In our derivation of dynamical friction, we neglected the boundary's curvature. It is straightforward 
that for a perturber on a normal trajectory that moves away from a convex boundary, the friction force 
will be further reduced with respect to that given by (|16ll7p as the wake behind the perturber is truncated 
further (as in Figure 1). 

The current analysis may also serve as a guide to examine the effect of the ambient medium's boundary 
in collisionlcss systems. As mentioned earlier, the derivation of dynamical frictio n from the velocity impuls e 
caused by a test star on the perturber is centered physically around the latter ( Binnev fc Tremamel 2008 ) . 



When steady state is considered, the large distance cutoff is applied in the perturber-centered frame which 
is equivalent to assuming a comoving symmetric boundary. Although they differ in magnitude, the friction 
forces in a gaseous m edium and a collisionless star system sh are a strong similarity in the way they depend 



on the Mach number (jSanchez-Salcedo fc Brandenburgll200ll ). Such similarity suggests that the conclusions 



we found here will hold for collisionless systems. 
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A. Gas flow near the medium's boundary 



In this Appendix, we characterize the flow in the presence of a reflecting planar boundary through the 
density enhancement and the velocity field. Prior to discussing the case of a medium with a boundary, we first 
determine the density enhancement and velocity field in an infinite medium. Much work has been devoted 
to the s tudy of a similar configuration that of the Bondi-Hoyle-Lyttleton accretion problem (for a review see 
Edgarl (|2004l) ). However most analytical studies of the flow have always dealt only with the ballistic flow 
for supersonic perturbers therefore excluding the effect of pressure in the medium. As explained below the 
ballistic flow occurs outside the wake region and is therefore irrelevant to dynamical friction because it is not 
accompanied by a density enhancement. In the context of Bondi-Hoyle-Lyttleton accretion with pressure, 
the analytical stud ies of the flow have characterized the veloc ity only in the vicinity of the singular and 



asymptotic points (jWolfson 



1977 



Yabushita 



1978: 



Horedt 



200 



0). 



The medium's response to the excitation of a point-like pcrturbcr of gravitational potential (j) p 
that is set in motion at time t = is described by the standard fluid equations: 



dtp + V • (po) = 0, 



d t v + (v ■ V)v = -- Vp - J^(t)V(f> p , 
P 



-GM/r 



(Al) 



where v(x,t) and p(x,t) are the perturbed velocity and density fields, p is the pressure, and Jf(t) is 
the Heaviside function. Prior to the perturber's motion, the gaseous medium is assumed to be at rest 
[i.e.v(x,t) = 0], homogeneous of density po, pressure po and isothermal with a sound speed c = (po/po) 1 ^ 2 - 
The linearization of medium's equation in the vicinity of the rest state yields: 



d t g + cV-v = 0, -dtu + Vg = -^P-Vc 



■pi 



(A2) 



where we used p = po(l + g) and v = cu. The combined equations in turn yield the forced sound wave 
equation satisfied by the density perturbation g as: 



1 



JT(t), 



V 2 g-^d 2 g = -^f>V 2 ^ p . 



(A3) 



The solution of this equation is known in electromagnetism theory as the Lienard-Wiechert potential and 
represents the electric potential of a moving electric charge distribution. It is written in integral form using 
the retarded Green function as: 



g(x,t) 
g(x,t) 



1 



+ 00 



dVdu 



6[u-t + \x- x'\/c] V 2 cj> p [x' - £(u)] 



AlTC 2 ./-oo \x — x 

f + °° S[u-t+\x- £{u)\/c] Jf(u) 



du. 



(A4) 



This preprint was prepared with the AAS IATjTjX macros v5.2. 
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The momentum equation (|A2[) shows that the velocity, v, can be written in terms of a potential function tp 
as v — given that the medium is at rest prior to the pcrturbcr's motion. The potential tp satisfies the 
equation: 



d t tp(x,t) 



1 



c 2 g(x,t) +(/> p (x,t)], 



which upon integration yields: 

tp(x,t) = - 

i> e (x,t) = - 

ip e {x,t) = 



GM 



c g(x, w) + 4> p (x, w)] dw 
o c 



(A5) 



(A6) 



* S 
du dw— 

u— — oo J w—0 



u - w + (R 2 + ((u) 2 ) 1/2 /c 



[R 2 + ((u) 2 } 1/2 
du 



0<«+(fl 2 +C 2 ) 1/2 /c<t, u>0 



y log ( z - Vu + [i? 2 + (z - Fit) 21 



dw 



o [i? 2 + CH 2 ] 



1/2 



los 



'z-Vf + \R 2 + {z-Vt) 2 ] 1/2 ' 



z + [R 2 + z 2 ] 



1/2 



(A7) 
(A8) 



wher e = z— £(«) and is V'baii the ballistic potential associated with the pressureless flow (jBisnovatvi-Kogan et al 
19791 ). To determine the boundaries of the density perturbation's contribution tp e to the velocity potential 



tp, we find the roots of cu + \J R 2 + £ 2 < ct subject to the conditions < u < t. These roots are given as: 

ct - J(z ± [(ct - Jtzf - (1 - M 2 )(ct 2 - r 2 )] 1/2 



u± 



c(l - Jt 2 ) 



(A9) 



For subsonic motion and outside the sonic wavefront, r = ct, there are no roots as there is no density 
perturbation. This implies that tp B = and tp — V'bali- Inside the sonic wavefront, r < ct, w_ is the only 
possible root so that both tp g and ipb&u contribute to tp and yield: 



tp(x,t) =-ylog 



-Vt+ [R 2 + (z- Vt) 2 } 



1/2 



V'sonic i) • 



(AlO) 



z - Vu_ + [R 2 + (z- Vu-) 2 ] 1 ' 2 , 

For supersonic motion, there are three domains: (i) outside the sonic wavefront and the Mach cone: ct < r, 
ct < (z + Ry/M 2 — V)j d% and ct/ 1 M < z < Met, (ii) inside the sonic wavefront: r < ct, (hi) outside the 
sonic wavefront but inside the Mach cone: (z + i? vC*^ 5 — 1 ) / j& < ct < r and ct/M < z < Met. Outside 
the sonic wavefront and the Mach cone, no roots are possible, tp e = and tp = Vwi- Inside the sonic 
wavefront, only «_ is possible and tp = tp sonic . Inside the Mach cone where (z + R\J Ji 2 — \)jM < ct < r 
and ct/M < z < Met, the two roots u± are possible and the velocity potential is given as: 



tp(x,t) = -prlog 



z-Vu++ [R 2 + (z - Vu+) 2 ] 1/2 



z-Vt+ [R 2 + {z - Vt) 2 ] 1/2 



Vu- + [R 2 + (z- Vu-)' 



,1/2 



z + [R 2 + z 2 } 



1/2 



lpMa.ch(x,t). 



(All) 



The calculated flow (|A8IA10IAlip is shown in Figure 8. The ballistic flow outside the wake region, V'baii, 
is finite despite the absence of a density enhancement. This results from the infinite interaction range of 
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the gravitational potential <p. This part of the flow however does not give rise to dynamical friction. The 
wavefront as a whole propagates as a shock in the gaseous medium. Inside the sonic wavefront, the velocity 
vanishes along the perturber's trajectory. For supersonic motion, whereas the velocity potential is continuous 
at the sonic wavefront r = ct, the velocity is discontinuous reflecting the discontinuity of the density g (|12l) . 
The charateristics of the flow given by ^ S O nic , and ^Mach agree with those obta ined from the numerical 



simulation of the nonlinear fluid equations f (jSanchez-Salcedo fc Brandenburall999| ). Figure 1). The velocity 



potentials we derive d does not agree with the expressions of the velocity field quoted without derivation by 
(l.lust fc Kege]ll990h . 



We now turn to the problem of a perturber on a trajectory £(£) in a medium with a infinite planar 
boundary defined as x ■ n = where n is its constant normal vector. The coordinate frame is chosen so 
that the boundary contains the origin. The medium fills the space where x ■ n > 0. For the outgoing wake, 
the density and velocity perturbations are truncated geometrically by the boundary as shown in Figure 1 
(bottom right panel). To obtain the density and velocity perturbations of the reflected wake, we use the 
property that normal fluid motion vanishes on the boundary. This may be written as v(x,t) ■ n = for 
x ■ n = 0. The equations of the density perturbation and the velocity potential to be solved for x ■ n > 
are given by: 

1 1 GM /"* 

V 2 gi g d\gx = - — V 2 p , -01 = / (cQi + <f> P )&T, 

c c c J 

and V^i ■ n = at x ■ n = 0. (A12) 



This problem can be solved by the method of images used in electromagnetism theory. The system that 
consists of a perturber and a boundary is equivalent to that of two identical perturbers placed on each side 
of the boundary and whose trajectories are mirror images of one another with respect to the boundary. To 
see this we note first that we seek the solution of the forced wave equation (|A12[) in the region x ■ n > 0. The 
second perturber, whose trajectory $, r (i) = — 2[£(i) -n]n, is located in the region x n < 0. Consequently, 
its corresponding density perturbation satisfies the force free wave equation in the region x ■ n > and could 
simply be added to the wave equation. Second, the trajectories of the two perturbers, £(t) and £ r {t), being 
mirror images of one another with respect to the boundary, the normal components of the respective velocity 
fields are antisymmetric and sum up to zero at the boundary. For a perturber moving inside the medium 
normal to its boundary initially at a distance d from it, we choose n = e z . We then have £(t) = (Vt + d)e z , 
£ r (t) = — (Vt + d)e z and the density perturbation and velocity potential are given as: 

qi = Jf(z)[g(R,z-d,V,t) + g(R,z + d,-V,t)], (A13) 
0i = Jf(z)[ip(R,z-d,V,t)+ip(R,z + d,-V,t)], (A14) 

where g and ip are given respectively by equations (fT2]) and (|A8l IA10I IA11|) . The first term in each equation 
corresponds to the geometrically truncated outgoing wake. The second term corresponds to the reflected 
wake. Figure 9 shows the surfaces of equal density and the velocity field for subsonic and supersonic 
perturbers. The reflected wake overtakes the subsonic perturber after a time t = 2d/(l — ./#). For a 
perturber moving inside the medium parallel to its boundary, we choose n = e x and £(t) = de x + Vte z . The 
mirror trajectory £ r (t) = —de x + Vte z and the density perturbation and velocity potential are given as: 

Pi = Jt?(x) \g{[{x-d) 2 + y 2 ]^,z 1 V 1 t S j + g ([{x + d) 2 + y 2 ]K z , V, , (A15) 

= Jf(x) [ijj([(x-d) 2 +y 2 ]?,z,V,t) +^([(x + d) 2 +y 2 ]Kz,V,t^ . (A16) 

Figure 10 shows the corresponding surfaces and the velocity field for subsonic and supersonic perturbers 
in the plane y — 0. Unlike motion normal to the boundary, the truncation of the wake is substantial. For 
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t 3> d/c the outgoing wake is reduced to approximately half its size in an infinite medium while the other half 
becomes the reflected wake. The presence of the two wakes on the same side of the pcrturber's trajectory 
gives rise to a force normal to the direction of motion pushing the perturber away from the boundary. 



0.5 1 1.5 -0.5 0.5 1 1.5 

z z 

Fig. 1. — Surfaces of equal density (|12[) for a subsonic perturber = 0.3) in a medium with a comoving 
symmetric boundary (left column) and a medium with a boundary unrelated to the perturber (right column) . 
Distances are scaled to r max . At t = r max /2c, the sonic wavefront (bold line) launched at t = has not 
reached the boundary (top row) and the forces in the two media are identical and given by (|5I9[) . In the 
steady state regime, t > 2r max /c (bottom row), the density perturbation is symmetric with respect to the 
perturber for the medium with the comoving symmetric boundary and the perturber is not subject to friction 
(bottom left panel) feciuations l5lll[) . The boundary of the stationary medium (bottom right panel) truncates 
the density perturbation asymmetrically and results in finite force on the subsonic perturber. 
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Fig. 2. — Friction force from the outgoing wake for motion normal to the boundary. The force scaled by 
2ir(GM) 2 po/c 2 is given by 7/./# 2 (|16I17|) ; time is scaled by d/c. The force is shown as a function of time 
for subsonic perturbers, M = to 0.9 (top left), supersonic perturbers, = 1.2 to 3.2, (top right) and 
for subsonic motion towards the boundary ^# = —0.9 to —0.1 (bottom left). The limit force as t ^ d/c is 
shown a function of the Mach number (bottom right) along the force experienced in an infinite medium Q 
(dotted line) for r Q = 10~ 2 d and ct/d = 10 2 , 10 3 and 10 4 . 
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Fig. 3. — Friction force from the reflected wake for motion normal to the boundary. The plots' parameters 
are identical to those in Figure 2. 
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Fig. 4. — Total friction force from outgoing and reflected wakes for motion normal to the boundary. The 
plots' parameters are identical to those in Figure 2. 
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Fig. 5. — Friction force fr<§m the Outgoing w3ke for motion parallel to the boundary. 4 The v| foifce scaled 
by 2tt(GM) 2 po/c 2 . The component along the trajectory (left column) is given by F 2 (ro/|l + j#\ < r < 
d) + F z (d < r < ct) where the first term is obtained from ([9j and the second term from (|3Tj) . The component 
normal to the trajectory (right column) is F x (d < r < ct) given by ([30)) ; time is scaled by d/c. The force is 
shown as a function of time for subsonic perturbers, ^ = to 0.9 (top row), supersonic perturbers, ./# = 1.2 
to 3.2, (middle row). The limit force as t 3> d/c is shown a function of the Mach number (bottom row) for 
r = 10~ 2 d and ct/d = 10 2 , 10 3 and 10 4 . 
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Fig. 6. — Friction force from the reflected wake for motion parallel to the boundary. The plots' parameters 
are identical to those in Figure 5. 
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Fig. 7. — 'Friction 1 force from the Sutgoing^ and" reflected w'akes for 1 motion 2 parallel 3 to the boufftiafy. The 
plots' parameters are identical to those in Figure 5. In the bottom left panel, the force along the trajectory 
for an infinite medium (|9|) is shown for comparison (dashed line) . In the bottom right panel, both F x (solid 
line) and \F Z \ (dashed line) are shown as functions of the Mach number. 
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Fig. 8. — Gas flow induced by a pcrturber with uniform rectilinear motion is an infinite gaseous medium. 
Distances are scaled to V/c and velocities in each panel are magnified by a constant factor to show their 
relative magnitudes. Top left panel: the sonic wave front of a subsonic pcrturber with ^# = 0.5 travels as a 
shock in the ambient medium while the outside pressurelcss ballistic flow points toward R = and z = Vt/2 
because of the symmetry of the corresponding potential. Top right panel: a zoom inside the sonic wave 
front around the pcrturber with ^ = 0.5. The velocity field is symmetric with respect to the trajectory 
and vanishes along it. Bottom left panel: the velocity field of a perturber that travel near the sound speed 
(.-# = 0.99) is strongly peaked at the perturber's position. Bottom right panel: For a supersonic perturber 
(jtft = 2), the velocity is maximal at the interface with the Mach cone. 
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Fig. 9. — Density enhancement (left column) and velocity field (right column^ for motion normal to the 
boundary at z = 0. Distances are scaled to d and time to d/c. Top row: A subsonic perturber with ./# = 0.5 
at t = 2; the reflected wake is behind the perturber. The solid lines denote the wavefront. Middle row: 
The reflected wake has overtaken the subsonic perturber (t = 5). Bottom row: A supersonic perturber with 
J( = 2. 
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Fig. 10. — Density enhancement (left column) and velocity field (right column) for motion parallel to the 
boundary at x = in the plane y — 0. Distances are scaled to d and time to d/c. Top row: A subsonic 
perturber with ^ = 0.5 at t = 2; the reflected wake is behind the perturber. The solid lines denote the 
wavefront. Middle row: The reflected wake has overtaken the subsonic perturber (t = 5). Bottom row: A 
supersonic perturber with ^# = 2. 



